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Abstract 

In axion quintessence, the cosmological era with an energy contrast 
in dark energy 0.1 < i^DE < 0.9 may represent a significant fraction 
of the universe's hfetime if the minimum of the axion potential is 
negative (unstable axion quintessence), thus resolving the cosmic co- 
incidence problem, as pointed out by Kallosh, Linde, Prokushkin, and 
Shmakova. Further details of the evolution of the quintessence field, 
the scale factor of the universe, and the Hubble parameter are pre- 
sented here, focussing on models with 0,ue,o = 0.72 and recent dark 
energy average equation of state parameter —l<wo< —0.85. For 
these parameter values, the contracting universe enters a late time 
era of kination, the negative Hubble parameter acting like a negative 
friction term in the Klein-Gordon equation, and the axion field makes 
many transits of — but never remains in — its vacuum state. 

Robust, scaled cosmological equations are derived for simulating 
the evolution of the scalar field, the scale factor, and the Hubble pa- 
rameter during a contracting phase of the universe. These equations 
allow the simulations presented here to proceed much closer to the 
singularity at the end of the collapsing universe than any previous 
simulations. 



1 Introduction 



Typically in quintessence theories with an asymptotically vanishing effective 
cosmological constant, the energy contrast in dark energy Qde rises from 



near zero for redshifts z > 5 to near one for z < —0.5, mimicking a true 
cosmological constant. At late times the quintessence field may begin to 
oscillate about its minimum, behaving like nonrelativistic matter, or the 
quintessence field may evolve toward infinity — in both cases with vanishing 
vacuum energy. In such theories, there is a period between roughly 3.5 Gyr 
and 20 Gyr after the big bang when 0.1 < ^Ide 0.9. However if the universe 
continues to expand forever, or even if positive curvature begins to dominate 
at late times (after the quintessence field has evolved to its minimum) and 
the universe enters a contracting stage, this period when the energy densities 
of dark energy and matter are comparable is a small or vanishing fraction 
of the total lifetime of the universe. This is called the cosmic coincidence 
problem. 

However in axion quintessence (as in other unstable de Sitter quintessence 
models), the cosmological era with 0.1 < Qde < 0.9 may represent a signifi- 
cant fraction of the universe's lifetime if the minimum of the axion potential 
is negative (unstable [de Sitter] axion quintessence), thus resolving [1] the 
cosmic coincidence problem. (Negative p\ or p^E and the fate of the uni- 
verse are discussed in Refs. [2] and [T] plus references therein.) 

Vilenkin [3] summarizes the currently predominant view that the dark 
energy is a cosmological constant, and points out the major weaknesses of 
most quintessence models: (i) most quintessence models assume that Pde ^ 
as the quintessence field evolves toward its minimum (which may be at 
\(j)\ oo), conflicting with the expectation that pde ought to evolve to a 
nonzero vacuum energy density; (ii) most quintessence models do not solve 
the coincidence problem that Pde,o = Pa ~ Pmo (the subscript "0" will denote 
present values); and (iii) since the present dark energy average equation of 
state parameter wq ~ —1, perhaps Wq = —1 and quintessence models are 
irrelevant. 

The unstable axion quintessence potential V{(f) = Acos{(f), where ip = 
(f)/Mp and the Planck mass Mp = I/VSttG = 2.4 x 10^^ GeV, addresses 
all of these issues, since the facts that the minimum of the potential is at 
—A ^ — Pa, Pde,q = Pa, and Wq ^ —1 but —1 are interrelated aspects of 
the model, and occur for an appreciable range of initial values for 0. 

For ^^(v?) = v4cos((y9), the initial value of the scalar field need only satisfy 
< v^i/vr < 0.23 to produce a universe like ours |1] (due to symmetry, we can 
restrict our attention to < < tt). Thus there is a significant 23% range 
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of the possible initial values which will produce a universe like ours|l] For 
these initial values, the contracting universe enters a late time era of kination 
(during which the scalar field kinetic energy dominates over all other forms of 
energy), the negative Hubble parameter acting like a negative friction term in 
the Klein-Gordon equation, and the axion field makes many transits of — but 
never remains in — its vacuum state|§ 

In Section 2, the basic cosmological equations are presented for the evo- 
lution of the scalar field, the scale factor, and the Hubble parameter, and 
cast in the form of a scaled, dimensionless system of first-order equations in 
the conformal time, appropriate for a contracting (or expanding) universe. 
These equations allow the simulations presented in Section 3 (see Figs. EHH]) 
to proceed much closer to the singularity at the end of the collapsing universe 
than any other simulations presented in the literature, and provide the basis 
for a more detailed analysis of the last stages of the collapsing universe than 
has appeared before. 



2 Cosmological Equations 

In the quintessence/cold dark matter (QCDM) model, the total energy den- 
sity p = Pm + Pr + P(j>, where pm, Pr, and are the energy densities in 
(nonrelativistic) matter, radiation, and the axion quintessence scalar field (j), 
respectively. Ratios of energy densities to the critical energy density pc for a 
fiat universe will be denoted by flm = Pm/Pc = Pr/Pc and = p^/pc, 
while ratios of present energy densities pmo, Pro, and p^o to the present crit- 
ical energy density pco will be denoted by Qmo, ^ro, and Q^o, respectively. 
flDE will denote 

_j ACDM 

~\n^ QCDM if < -1/3. 

Using WMAP5 [S] central values, we will set Qde,o = 0.72, •0 — 0.5 X 
10-^ Qmo = I - ^DE,o - ^ro ~ 0.28, and p^^^ = 2.5 x IQ-^ eV, with the 
present time to = 13.73 Gyr after the big bang for ACDM. 

^Qualitatively similar results to those presented here are obtained for V{(p) = A cos{Xlp) 
for A = 0(1). 

^The coupling of the quintessence field to other particles must be very small, and will 
for the most part be neglected in this investigation. 
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The homogeneous scalar field obeys the Klein-Gordon equation 



The evolution of the universe is described by the Friedmann equations for 
the Hubble parameter H = a/a and the scale factor a{t) 



p k 

1 



6M 



(P + 3P) (4) 



where the energy density p = + Pm + Pr and the pressure P = P(j, + Pm + Pi 
with Pm = 0, Pr- = Pr/3, and 



r-i 



P0 = ^0' + n0), P<t> = \^''-V{<i>). (5) 

The curvature signature k = +1, 0, —1 for a closed, fiat, or open geometry. 
Eq. (jlj) shows that P < — p/3 for an accelerating universe. 

The conservation of energy equation for matter, radiation, and the scalar 
field is 

p + 3H{p + P) = 0. (6) 

Equation ([6]) gives the evolution of pm and pr, and with Eq. ([5]) the Klein- 
Gordon equation ([2]) for the weakly coupled scalar field. The time rate of 
change of the Hubble parameter is given by 

Only two of Eqs. ([3]), @, IQ, and ([7]) are independent. We will assume a 
fiat universe after infiation and henceforth set k = 0. 

The logarithmic time variable (number of e-folds of the scale factor) is 
defined as r = ln(a/ao) = — ln(H-2;). Note that for de Sitter space r = H\t, 
where Hi = pa/(3M|,), and that i^A^ is a natural time variable for the 
era of A- matter domination (see e.g. Ref. P). We will make the simple 
approximations 

Pr = PrOe""^^, Pm = PmOe~^^ ■ (8) 
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The equation of state parameter for the scalar field is w,^ = P^j p^. 
Since r is a natural time variable for the era of A-matter domination, we 
define the recent average of as 

1 

wo = - w^dr. (9) 

T Jo 

We will take the upper limit of integration to correspond to z = 1.75. The 
SNe la observations [7J bound the recent average —1.1 < wq < —0.85 (95% 
CL). 

For numerical simulations, the cosmological equations should be put into 
a scaled, dimensionless form. Equations ([2]) and ([3]) can be cast [1] in the 
form of a system of two first-order equations in r plus a scaled version of H: 

= (10) 



= -W^ (11) 
= p (12) 

P=li^^ + V + Pm + Pr (13) 


where ip = ^/Mp, ^ = e^^ip/Ho, H = e^^H/Ho, V = e^^V/pco, % = 

e^'^V^/PcO, P = e^^P/PcO, Pm = e^'^PnJpcQ = r^mOe^, Pr = e'^'^Pr/PcO = ^rO- 

This scaling results in a set of equations that is numerically more robust, 
especially near and before the time of big-bang nucleosynthesis (BBN) — see 
Ref. [5, especially Fig. 1. 

For a contracting phase (in which H goes through zero), a different set 
of equations and a different scaling should be used. Here we will use the 
conformal time variable 

V = r ^dt (14) 

JO a 

where t = corresponds to the big bang. 

Equations (E]) and (I7j) can be cast in the form of a system of three first- 
order equations in rj: 

^ = -2H^Ij - (16) 

077 
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dH 1 / . „— \ 1 2 , Y7 ^mO-r n .-2t 



:^ = --{p + 3P)=--r + V- ^e-^ - n^oe-'^ (17) 
P = l^P' + V + p^ + p,., P=lv;2-F+^p, (18) 



where ip = (f)/Mp, ^/j = e^ip/Ho, H = e^H/Ho, V = e^^V/pco, = 
e^W^/pcQ, p = e^^p/pco, Pm = e^^pm/pco = ^moe"^, Pr = e^^pr/pco = 
e~^^f2ro- This scahng resuhs in a set of numerically more robust equations, 
especially near the turn-around time between expanding and contracting 
phases of the universe. 

Note that the conformal time t] is related to the logarithmic time r by 

dr — , , 

3 Simulations of Unstable Axion Quintessence 

The original axion quintessence potential V = A{1 + cos (99)) was based on 

= 1 supergravity [H [9], with = 3H^. As vr, the universe evolves 
to Minkowski space. 

The unstable de Sitter axion potential V = y4cos(v9) is based on M/string 
theory reduced to an effective = 1 supergravity theory [TU], with = 
—3Hj^ at the maximum of V. 

Both axion quintessence models are derivable (up to a constant) from 
string theory as axion monodromy [llj. 

The quintessence axion is a pseudo Nambu-Goldstone boson: at the per- 
turbative level the theory is shift symmetric under ip ^ (f + const with 
(fi = (f)/M. The shift symmetry is broken — before or during inflation — by 
nonperturbative instanton effects to a discrete symmetry ip ^ (p + 27r, gener- 
ating a potential V{(f) = A{C + cos((^)). In these theories, quantum correc- 
tions to the classical axion potential are suppressed. For quintessence (or for 
natural inflation [12]), M ~ Mp; we will take M = Mp. C = and C = 1 
are the most interesting unstable axion and original axion cases, respectively. 

For the unstable axion quintessence computations (with an expanding 
and contracting universe), we will use Eqs. f|T5|) - f|T7|) and f|T9|) with initial 
conditions ipi and 0j = specified at matter-radiation equality Zmr = 3280, 
which corresponds to 

= /o^ /T^^ • (20) 
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The constant A in the potential is adjusted so that fi^o = 0.72. This involves 
the usual single fine tuning. Note that the same anthropic arguments that 
limit the magnitude of a present-day cosmological constant also limit A < 
lOOpcO; so the tuning of A is no worse than the tuning of a cosmological 
constant. 





A/pco 


'"^0 


to 


^0.1 


^0.9 










0.05 


0.73 


-0.998 


13.72 


3.6 


20.0 


63.2 


72.7 


0.23 


12.0 


0.10 


0.78 


-0.99 


13.70 


3.5 


20.3 


47.6 


56.8 


0.30 


5.0 


0.15 


0.88 


-0.97 


13.64 


3.5 


21.2 


37.6 


46.2 


0.38 


3.0 


0.20 


1.09 


-0.93 


13.49 


3.3 


21.2* 


29.2 


37.0 


0.48 


2.0 


0.23 


1.41 


-0.87 


13.25 


3.0 


16.8* 


23.9 


30.7 


0.45 


1.6 



Table 1: Parameters for the potential V = Acos{lp). to is the current age 
of the universe with to = 13.73 Gyr in the ACDM model, 0.1 < ^Ide < 0.9 
for to.i < t < ^0.9, the "coincidence" time interval Ate = ^o.g — ^o.i, t* is the 
turn-around time, t/ is the time of the big crunch, and a* = a{t^,). All times 
are in Gyr. *For (pi/n = 0.20 (0.23), ilnE < 0.85 (0.77) and in these cases 
^0.9 = to.85 and to.77, respectively. 

Results for the unstable axion potential are presented in Tabled] for var- 
ious ipi and for y^j/vr = 0.1 in Figs. HHHl (As — )■ 0, classically tf ^ oo, but 
quantum effects destabilize ~ so that the maximum t/ ~ 100 to [I]-) For 
the values in the Table, as ipi increases, (0j)| also increases and starts to 
move earlier, leading to a decrease in to.i, ^O; t*, and tf, and correspondingly 
to an increase in wq away from —1. Note that for (fi/ir = 0.2, the coincidence 
time ratio approaches 50%. 

The QCDM universe mimics the ACDM model (see Fig. [T| for clarity, 
only the beginnings of the contracting stage are shown in this figure) until 
about z = —0.5, after which the QCDM universe begins to decelerate and 
ultimately to rapidly contract to a big crunch (Figs. [SJandE]). 

During the contracting stage, H < acts as a negative friction in 

4> + 3H<i) + V^ = 

amplifying the axion field and its kinetic energy to bring about a late stage 
kination era during which the scalar field kinetic energy dominates over all 
other forms of energy. 
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The quintessence axion is an ultra-light scalar field with ~ H\, so 
"sits and waits" during the early evolution of the universe, and only starts 
to move when H"^ ~ (Figs. [Hand [5]). In this way it is easy to satisfy the 
BBN {z ~ 10^-10^^), cosmic microwave background (CMB) {z ~ 10^-10^), 
and large scale structure (LSS) {z ~ 10-10'*) bounds on VLde ~ 0.1, as in 
Fig. [U An ultra-light scalar field also refiects the observational evidence that 
the universe has only recently become dominated by dark energy. 

In Fig. [21 the Hubble parameter goes through zero at the turn-around 
time between an expanding and contracting universe. At the beginning of 
the contacting stage, has yet to reach the minimum of the potential energy 
(see Fig. E]), and thus the negative Hubble parameter amplifies the kinetic 
energy of the scalar field, bringing about an era of kination with = 1, as 
seen in Figs. IH O O and [HI Also note that Figs. [2l [31 and [5] indicate that H, 
0, and are approaching a singularity near tf. In fact, in an era of kination 



during contraction during which H = — (j) /(-\/6Mp), = ^Ij'ij {tf —t) 
and (f = -J2/3\n {tj - t), while a ~ (t/ - t)^^^ (see Fig. E>- 



Figure [6] shows that as increases without bound, the potential energy 
V{(f)), since a periodic function of 0, oscillates more and more rapidly. De- 
pending on the strength of the coupling between the quintessence axion and 
other particles, may decay and populate the universe with additional radia- 
tion and matter. (The monotonically increasing field in a periodic potential 
can be interpreted as an oscillating field.) At the very least, there should be 
gravitational production of particles by during contraction. 

Figures [3 and [8] follow the contracting stage further, illustrating that 
although — ±oo twice just after t^:, the universe ultimately enters a 
stage of kination in which W(i) = 1 near tf. Note that ~ —1 until well 
after to. 

After matter-quintessence equality at tm(f, ~ 0.7 to, the scalar field energy 
density always dominates over the matter and radiation energy densities. 
There is a period from t ^ 2.5 to until 3.8 to when the scalar field potential 
energy is comparable to its kinetic energy, and then the kinetic energy (which 
scales as 1/a^) predominates during the rapid contraction to a big crunch (see 
Fig.M)- 
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4 Conclusion 



As the universe contracts, density inhomogeneities are amplified and pre- 
sumably black holes are formed, similarly to the contracting stage of the 
ekpyrotic universe [131 E] with w = 1. Depending on the strength of the 
coupling (which we have neglected here) between the quintessence axion and 
other particles, (j) may decay and produce radiation and matter. There should 
at least be gravitational production of particles by (p during contraction. As 
the universe reheats during contraction, broken symmetries are restored. It 
is possible that inflating patches may be generated, spawning new universes 
from the old. These issues are currently under investigation. 

In summary, the unstable axion quintessence potential resolves the issues 
concerning quintessence raised by Vilenkin: the minimum of the potential is 
not at zero, but at a negative value ~ — Pa, Pde,o = Pa with only a single fine 
tuning (in the anthropic range), and Wq naturally satisfies —l<Wo< —0.87, 
for an appreciable 23% range of possible initial values for the quintessence 
field. And for a universe like ours, the coincidence time when the energy 
densities of dark energy and matter are comparable varies (as long as (pi/n 
is not too small — say, > 0.05) from 25%-50% of the lifetime of the universe. 
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Figure 1: Q vs. logiQ(l + 2) for the potential V = Acos{ip), (pi/n = 0.1 (solid) 
vs. ACDM (dotted). The light yellow rectangles are the bounds on Qde from 
LSS, CMB, and BBN. 
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Figure 8: Overview of vs. t/to. The dot indicates 




Figure 9: Scalar field kinetic energy density 0^/2 and potential energy density 
V{(j)) vs. t/tQ. The dot indicates t*. 
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